Mon. Not. R. Astron. Soc. 0Q0.mfl4l(2Q09) Printed 4 March 2010 (MN KTeX style file v2.2) 



New 9.9-GHz methanol masers 



o 



6 



00 

m 

p 

rn 
O 

o 



X 



M. A. Voronkov^'^*, J.L. Caswell\ S.P. Ellingsen^ A.M. Sobolev^ 

^Australia Telescope National Facility, CSIRO Astronomy and Space Science, PO Box 76, Epping, NSW 1710, Australia 
^ Astro Space Centre, Profsouznaya st. 84/32, 1 17997 Moscow, Russia 

^School of Mathematics and Physics, University of Tasmania, GPO Box 252-37, Hobart, Tasmania 7000, Australia 
"^Ural State University, Lenin ave. 51, 620083 Ekaterinburg, Russia 



ABSTRACT 

The Australia Telescope Compact Array (ATCA) has been used to make the first extensive 
search for the class I methanol masers at 9.9 GHz. In total, 48 regions of high-mass star for- 
mation were observed. In addition to masers in W33-Met (G12. 80-0.19) and G343. 12-0.06 
(IRAS 16547—4247) which have akeady been reported in the literature, two new 9.9-GHz 
masers have been found towards G331. 13—0.24 and G19.61— 0.23. We have determined ab- 
solute positions (accurate to roughly a second of arc) for all the detected masers and suggest 
that some class I masers may be associated with shocks driven into molecular clouds by ex- 
panding Hll regions. Our observations also imply that the evolutionary stage of a high-mass 
star forming region when the class I masers are present can outlast the stage when the class II 
masers at 6.7-GHz are detectable, and overlaps significantly with the stage when OH masers 
are active. 

Key words: masers - ISM: molecules - ISM: Hll regions 



1 INTRODUCTION 



Methanol masers are well established probes of high-mass star 
formation. They a re divided into two categories first defined by 

iBatrla et al. I l ll987h : class II masers (of which the 6.7 GHz is the 
best known and usually st rongest) are ass ociated with millimetre 
and infrared sources (e.g.. lHilleta"i]|2005t) and reside in the close 
environment of high-mass young stellar objects (YSOs). Class I 
masers (of which the 44 and 95 GHz are commonly observed) are 
' often found apart from the strong continuum and infrared sources 
and can be separated by u p to a p arsec from the YSO re spon- 
sible for excitation (e.g ., IXurtz et al.ll2004l : IVoronkov et al.ll2006l : 

ICyganowski et alj2009() . 

Theoretical calculations are able to explain this empirical 
classification and strongly suggest that the pumping process of 
class I masers is dominated by collisions with the molecular hydro- 
gen ( see, e.g.. lLees|[l973l : ISobolev & Strelnitskiilll983l : IVoronkovl 

1 19991) in contrast to class II masers, which are pumped by radia- 
tive excitation (see, e.g., IWilson et alj|l985l : ISobolev & Deguchil 

1 1994 [ Sutton et al. The t wo pumping mechanisms were 



shown to be competitive (see, e.g ., Mentenll99lMCragg et al.l 19921 : 



IVoronkovlll999l : IVoronkov et al.ll2005bir For example, strong radi- 
ation from a nearby infrared source quen ches class I masers an d 
increases the strength of class II masers dVoronkov etai]|2005bh . 
Although it has been demonstrated by theoretical calculations that 
a weak class II maser at 6.7 GHz c an coexist with a brigh t class I 
emission under special conditions jVoronkov et alj|2005bt) , bright 
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masers of different classes residing in the same volume of gas are 
widely accepted as mutually exclusive. However, on larger scales, 
they are often observed to coexist in the same star forming region 
within less than a parsec of each other. 

Class I methanol masers are relatively poorly studied. The 
common consensus is that the majority of class I masers trace inter- 
face regions between outflows and molecular gas, although direct 
observational evidence of t his has been obtained for a l i mited num- 
ber o f sources only (e.g., Plambeck & Mentenl 1 19901 : iKurtz et al] 
l2004l : IVoronkov et alj |2006|) . An alternative scenario, which in- 
volves cloud-cloud collisions, may be realised in some sources 
(Sobolev 1992; Mehringer & Menten 1996; Salii et al. 2002). The 
common ingredient of these two scenarios is the presence of 
shocks. There is observational evidence that the methanol abun- 
dance is significantly increased in the sh ock processed regions (e.g. 
lGibb&Davi3ll998l : ISutton et al.ll2004h . The gas in such regions 
is heated and compressed increasing the frequency of collisions 
with molecula r hydrogen and, th e refore, prov i ding a more efficient 
pumping te.g. lSutton et al.ll2004) . IChen et all (12009) demonstrated 
statistically the association of class I masers with the shocks traced 
by extended features showing a prominent excess of the 4.5-^m 
emission in the images obtained with the Spitzer Space Telescope's 
Infrared Arr ay Camera (IRAC) also k nown as extended green ob- 
jects tEGOs: ICvganowski et al.l2008l) . 

To date, there is observational evidence that more than 20 
methanol transitions can produce class I masers, more than half 
of which belo ng to the J2— Ji E series in the 25 —30 GHz fre- 
quency range jVoronkov! [T999I : iMuUer et al.ll2004) . The observa- 
tional properties of different class I masers are not the same. Some 
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masers, e.g. at 44 and 95 GHz, are known to be quite ubiquitous 
with a bout 100 sources k nown to date (see, e .g., |H aschick et aL, 
1990'; 'Slvsh et al.' 'l994'; IVal'tts et alj l200d : iKurtz et ali 120041 
Cyganowski et al. 2009). On the other hand, there are other tran- 
sitions with very few known maser sources. An example of such a 
transition showing rare maser activity is the 9_i— 8_2E methanol 
transition at 9.9 GHz. The rarity of these masers is likely due to the 
strong dependence of the maser brightness on the physical condi- 
tions, especially in the requirement of higher typical temp eratures 
and densities for these masers to form jSobolev et alj200a) . 

ISlvshetallj 1993 1) conducted the only search for the 9.9-GHz 
masers reported in the literature so far. However, they observed just 
11 targets which were largely well known regions of high-mass 
star-formation where methanol masers of either class had previ- 
ously been reported and discovered a single 9.9-GHz maser (W33- 
Met also known as G12.80— 0. 19). In this paper we report the re- 
sults of the search for the 9.9-GHz masers towards 48 independent 
positions. The majority of targets, 46 sources i n total, were know n 
class I masers at 44 and/or 95 GHz dra wn fr omlSlvsh et al.Ul994h . 
IVal'tts et alj ( l2000t) . lKurtz et al.l ^20041) . and lEUingserJ ( l2005h . The 
remaining two sources were included in order to cover all known 
star forming regions located south of declination of —20°, which 
have a per iodically variable class II methanol maser at 6.7-GHz ac- 
cording to lGoedhart et 'Zll l2004h . Little is known about the physics 
of periodic variability of such masers. Therefore, a detection of the 
9.9-GHz maser (which is very sensitive to the physical conditions 
and has a different pumping to the 6.7-GHz maser) in the same 
source might shed light on this enigmatic phenomenon. 



2 OBSERVATIONS 

Observations were made with the Australia Telescope Compact Ar- 
ray (ATCA) in 7 time allocations scheduled from December 2005 
to July 2008 (the archive project code is C1466). The details of 
these observing sessions along with the array configurations used 
and the range of baseline lengths are summarised in the first five 
columns of Table[T] Each source was observed in several (typically 
6) cuts of 10 minutes each to get a good hour angle coverage. The 
data reduction was performed using the MIRIAD package (12 June 
2007 release) following standard procedures (with the exception of 
bandpass calibration). The bandpass calibration was achieved by 
fitting a low-order polynomial to the spectra of the bandpass cali- 
brator (1921-293, except for the observations on December 17 and 
June 2, when 0537-441 and 1253-055 were observed, respectively) 
using the UVLIN and MFCAL tasks. This approach for bandpass cal- 
ibration is usually superior and adds less noise than an independent 
solution for each spectral channel, if one tries to detect a weak nar- 
row line in the presence of a strong continuum, which is the case 
for many sources in our sample. 

The regular observations were preceded by a short test project 
(the archive project code is CX075) using Director's time in De- 
cember 2004 and January 2005. It aimed at testing the viability of 
observations at 9.9 GHz, which is outside the nominal range of the 
ATCA X-band receivers (8.0 to 9.2 GHz), but yielded the serendip- 
itous discovery of a new 9.9 GHz maser in G343. 12-0.06 (IRAS 
16547-4247 ). This source was the n investigated in detail in a sep- 
arate paper jVoronkov etal]|2006l) . The receiver performance was 
found to vary significantly from antenna to antenna at 9.9 GHz due 
to the proximity of the band edge. In particular, baselines including 
antenna CAOl had approximately 3 times higher noise than other 
baselines in most sessions. As the masers are unresolved sources 



for the ATCA, this inhomogeneity in the sensitivity of individual 
baselines is taken care of automatically by the imaging process. 

We adopted an astronomical measurement of the 
rest frequency for the 9-i — 8-2 E m ethanol transition, 
9936.201±0.001 MHz j Voronkov et"al] I2OO6I) . The uncertainty 
corresponds to 0.03 km s~^ of the radial velocity. The adopted 
rest frequen cy lies within the u ncertainty of the laboratory mea- 
surement of iMiiller et al. I (12004) equivalent to 0.12 km s ^. The 
correlator was configured with 1024 spectral channels across a 
4 MHz bandwidth, providing a spectral resolution of 0.12 km s^^. 
In addition, the spare correlator capacity was used to observe 
the 8.6 GHz continuum emission (with the 128-MHz bandwidth) 
simultaneously with the line measurement. Although two or- 
thogonal linear polarisations have been recorded, the correlator 
configuration used in the project did not allow us to calibrate 
the instrumental polarisation and hence to determine polarisation 
properties of the observed emission. During the data reduction, the 
orthogonal polarisations were converted to Stokes-I. Therefore, all 
flux densities presented in the paper are not affected by polarisation 
to the first order. 

All image cubes were constructed from the visibility data af- 
ter continuum subtraction. The natural weights were used provid- 
ing a typical synthesised beam size of around a few seconds of 
arc. Due to software limitations the search area was split into 9 
sub-cubes with a 50% overlap. Each sub-cube containing 512x512 
spatial pixels of 0.5x0.5 arcsec and 1024 spectral channels was 
then processed separately. This strategy covered a rectangular area 
of 8.3 arcmin in size. The half power width of the primary beam 
is 5.1 arcmin at this frequency. An automated spectral line finder 
(Voronkov, in prep.) based on the ATNF Spectral line Analysis 
Package {http://www.atnf.csiro.au/computing/software/asap ) was 
used to search for sources in the dirty spectral cubesQThe spectral 
line finder applied a simple statistical criterion of at least 3 con- 
secutive spectral channels above the threshold of 2x7 (with respect 
to the noise in the local vicinity of the tested channels) to claim 
a detection. The 3-dimensional search routine formed slices along 
the spectral axis for each spatial pixel, ran the 1 -dimensional al- 
gorithm and combined the results into sources. The sources which 
had less than 4 spatial pixels with detected emission (half power 
area of the main lobe of synthesised beam) were discarded. Due to 
statistical noise fluctuations, the algorithm produced of the order of 
10 spurious sources per cube, each of which was examined manu- 
ally. The same source detection software was used and thoroughly 
tested in the Parkes 6-GHz Multibeam survey (hereafter MMB sur- 
vey; Green et al. 2009). In addition to automated software, both 
single baseline and averaged uv-spectra (visibility amplitude ver- 
sus frequency) were inspected by eye to search for possible broad 
thermal components, which may not be imaged adequately with the 
uv-coverage provided by the sparse ATCA configurations. Due to 
the low dynamic range in the dirty image cubes the automated soft- 
ware reported a number of detections around the field of view for 
each real maser. In this case, a deconvolved cube was made using 
the CLEAN task and the source parameters were measured using the 
IMFIT task of MIRIAD. A self-calibration was performed and a new 
cube was constructed after the absolute position had been deter- 
mined. The spectra were extracted from such self-calibrated cubes 
using the IMSPEC task. An additional search for weaker masers 
was performed in the parts of the self-calibrated cube unaffected 



^ Strictly speaking, sidelobes are only caused by real sources. Therefore, 
no deconvolution is required unless there is a detection. 
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Table 1. Dates of observations, a summary of array configurations and decorrelation estimates. More details on the ATCA configurations are available on the 
web {http://www.narrabri.atnf.csim.au/ohsennng/configs.htinI i. Note, that antenna CAOl did not perform well at this frequency in 2005/2006. It is involved, 
in particular, in the shortest spacing of the 6C an'ay configuration. Therefore, the shortest baseline was not as sensitive as the other baselines in March 
observations. Antenna CA06 did not observe in 2008 sessions. 



UT Date 


Array 




Baseline length 


Decorrelation estimate 




of 


configuration 


min 


max 


max 


Source 


Separation 


Estimated 


observations 






no CA06 


with CA06 


and 


angular 


temporal 


factor 






(m) 


(m) 


(m) 


calibrator 


(deg) 


(min) 




15 Dec 2005 


6A 


337 


2923 


5939 


1059-63 (odd+even) 





23.8 (2) 


3.8 (7) 












1613-586 (odd-l-even) 





21 (5) 


5.7 (5) 


16 Dec 2005 










1613-586 (odd-l-even) 





18 (3) 


6.8 (4) 


17 Dec 2005 










1613-586 (odd-l-even) 





17 (3) 


1.88 (3) 


25 Mar 2006 


6C 


153 


2786 


6000 


1646-50 (odd-l-even) 





27.8 (1) 


1.179 (4) 












1710-269 -H 1730-130 


14.6 


20.1 (7) 


1.25 (1) 


26 Mar 2006 










1646-50 (odd-l-even) 





20.4 (2) 


1.14 (2) 












1710-269 (odd-l-even) 





21 (1) 


1.250 (9) 


2 Jun 2008 


EW352 


31 


352 


none 


1646-50 (odd-l-even) 





17 (2) 


1.004(5) 


20 Jul 2008 


H214 


92 


247 


none 


1646-50 (odd-l-even) 





9 (2) 


1.036 (2) 












1613-586 (odd-l-even) 





9 (2) 


1.065 (3) 












1646-50 + 1613-586 


9.3 


1.8(1) 


1.16 (2) 



by the main maser feature using the same approach as for the orig- 
inal cubes. However, we did not find any additional weak maser 
features. 

During the first two observing allocations of the regular ob- 
servations in December 2005 the atmospheric conditions were very 
unstable with thunderstorms present in the area causing an interrup- 
tion once on December 15. Observations through adverse weather, 
aggravated by the sparse array configuration, resulted in an in- 
creased level of decorrelation due to a poor phase stability. As 
this decorrelation directly affects the flux density limit for all non- 
detections, we have attempted to estimate it using the calibrator 
data. We used two different methods, depending on the data avail- 
able for a particular observing session. On March 25 two separate 
continuum calibrators 1710-269 and 1730-130 were observed us- 
ing the same frequency setup, which was chosen to achieve an ap- 
propriate velocity range for the target sources. The first source of 
this pair acted as an ordinary secondary calibrator allowing us to 
form a calibration solution. The second source was then imaged 
using this calibration. Another image of the second source was con- 
structed following a self-calibration. Then, the decorrelation factor 
was measured as the ratio of the two flux densities obtained from 
each image using the IMFIT task. The second method, suitable for 
observations with only one calibrator, is to split the whole sequence 
of calibrator scans into odd and even scans and then treat the odd 
ones as a calibrator and the even scans as a source. The decorrela- 
tion factor was then obtained from the images the same manner as 
for the first method. 

The right-hand four columns of Table [T] summarise the re- 
sults for each observing sessions, including the method and sources 
used, angular and mean temporal separation between the would be 
target and calibrator, and the estimated decorrelation factor. The 
uncertainties are expressed in the units of the least significant fig- 
ure. Table [T] shows that the decorrelation was up to 7 times worse 
during the first two observing sessions. It is worth noting that these 
estimates are likely to be larger than the actual decorrelation af- 
fecting the target sources because we effectively double the length 
of the calibrator cycle when splitting calibrator scans into odd and 
even scans. For some observing sessions more than one decorre- 
lation estimate is available (Table [TJ. The estimates are in a good 
agreement for the session on December 15, when the phase stability 



was bad, but are more than 3a apart for the March sessions, which 
had a good phase stability. This is likely to be due to systematic 
effects such as different declinations for the calibrators involved in 
the estimates (which makes the average elevations different), and 
an additional angular separation when a pair of sources was used 
instead of the odd-l-even method. We used the largest decorrelation 
estimate obtained for each observing session to scale the flux den- 
sity limits for non-detections. 

The last two sessions in 2008 (Table [TJ were completely de- 
voted to the 9.9 GHz maser discovered towards G33 1 . 1 3 — 0.24 dur- 
ing the session on December 16 and no new targets were observed. 
The goal of these repeated observation was to provide an indepen- 
dent measurement of the absolute position of the maser. We found 
that the morphology of G3 3 1.1 3— 0.24 was such that an indepen- 
dent estimate of a systematic uncertainty was required for an un- 
ambiguous interpretation. A typical (la) astrometric accuracy of 
ATCA (largely determined by systematic effects for masers) is usu- 
ally adopted to be around 0'.'5. However, the bad weather of the De- 
cember 16 session is likely to have further degraded the accuracy. 
Although two 2008 sessions used relatively compact array configu- 
rations and the synthesised beam was an order of magnitude larger, 
by observing in good weather and by using a different calibrator 
(1646-50 instead of 1613-586) for which a more accurate position 
is available, we have achieved a comparable formal uncertainty of 
the position measurement in all three sessions. We found that the 
absolute positions obtained in all three sessions are consistent, and 
we estimate the final accuracy (sigma for the mean) to be 1 arcsec- 
ond. 

3 RESULTS 

In addition to the new 9.9-GHz maser found in G343.12— 0.06 
during the pilot survey as a chance discov ery (results of the 
follow-up are reported in ' Voronkov et al]|20 06'), the regular survey 
yielded three other detections, listed in Table |2] These include the 
known 9.9 GHz m aser in W33-Met (G12.80-0.19) discovered by 
ISlvshetalJ ( fl993h , and two new detections. The first six columns of 
the table are self-explanatory and include the date of observation, 
the synthesised beam size and its position angle, the noise level 
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Table 2. Observation details for sources detected at 9.9 GHz. All uncertainties are given in the brackets and expressed in the units of the least significant figure. 
The uncertainty of the decorrelation factor was calculated from the uncertainties of the flux measurement. 



Source 


Date of 


Synthesised beam 


la Noise 


Inspected 


Measured 


Separation 


from calibrator 


name 


observation 


FWHM 


p.a. 


level 


velocity range 


decollation 


angular 


temporal 






(arcsec) 


(deg) 


(mJy beam^ ^ ) 


(km s-i) 


factor 


(deg) 


(min) 


G331. 13-0.24 


16 Dec 2005 


2.1x1.5 


15 


39 


-124, -27 


2.9 (3) 


7.0 


5.3 (1) 




2 Jun 2008 


18x12 


3 


35 


-131, -34 


1.04 (2) 


6.2 


4.0 (1) 




20 Jul 2008 


17x14 


53 


68 


-150, -53 


1.09 (3) 


6.2 


4 (1) 


W33-Met 


25 Mar 2006 


2.6x1.0 


14 


44 


-29 ,+53 


1.18 (4) 


11.0 


5.4 (1) 


G19.61-0.23 


25 Mar 2006 


3.7x1.0 


7 


43 


-28 ,+68 


1.60 (6) 


13.4 


21 (3) 



in the image representing the detection threshold, and the veloc- 
ity range searched for maser emission. The last three columns are 
related to the decorrelation estimates. For sources with a maser de- 
tection, self-caHbration effectively provides a measurement of the 
decorrelation factor. The factor given in Table |2] was calculated as 
the ratio of the flux density obtained using a self-calibrated data set 
of the maser to that prior to the self-calibration. The values are be- 
low the estimates obtained for a particular observing session using 
the calibrator data (see Table [TJ as they are expected to be, with 
the exception of the last source, G19.61— 0.23. For this source, the 
measured decorrelation is slightly higher than the estimate deduced 
from the calibrator data. This fact can probably be explained by a 
more northern location of the source than either of the calibrators 
used for decorrelation estimates and, as a result, a lower elevation 
on average throughout the observing run. Systematic effects such 
as fit errors due to the elongated synthesised beam along the dec- 
lination axis, which is produced by the East- West arrays for the 
sources located near the equator, could also contribute to this minor 
discrepancy. For the purpose of comparison, the last two columns 
of Table |2] give the angular and temporal separation of the maser 
source from the appropriate secondary calibrator. The uncertainties 
are given in round brackets and expressed in the units of the least 
significant figure. 

The absolute position, peak velocity and flux density are 
summarised in Table [3] for all detect ed masers. We have added 
G343. 12-0.06 ( I Voronkov etai]|2006 l) so that Table [3] contains a 
complete list of all four 9.9 GHz masers known to date. The un- 
certainty of the peak velocity is the spectral resolution given in the 
previous section (the rest frequency uncertainty is much smaller 
than the spectral resolution). The flux density uncertainties given 
in Table [3] are cumulative values including both the image noise 
given in Table |2] and the accuracy of the flux scale calibration, 
which is typically of the order of 3% at cm-wavelengths. The spec- 
tra of the 9.9-GHz detections are shown in Fig.[T] All masers have 
simple spectra in this transition containing just one component. In 
contrast, most sources have multiple components in the spectra of 
the most wi despread class I methanol maser transition at 44 GHz 
dSlvsh et al.lll994l : lKogan & Slvshlll998l) . This situation is similar 
to that observed in G343.12— 0.06, where the 9.9 GHz maser and 
other rare class I masers at 25 and 104 GHz are confined to one 
spot only, while the more widespread 95 and 84-GHz masers have 
been observed in a number of spots (Voronkov et al. 2006). Such a 
behavior confirms the hypothesis that the 9.9 GHz masers are pro- 
duced under a more restricted range of physical conditions than the 
widespread masers. 

Figure[2]shows positions of detected 9.9-GHz masers with re- 
spect to the 8.6-GHz continuum observed simultaneously. The po- 
sitions of the 6.7-GHz masers (see notes on individual sources for 
references) are also shown. The grayscale represents the Spitzer 



IRAC 4.5-/^m image s. For one source, G 331.13— 0.24, we used 
the contiuum data of IPhillips et alj ( Il998l) . Both the morphology 
of the Hll region and the flux den sity obtained i n our measure- 
ment ar e in a good a g reeme nt with IPhillips et alj ( Il998[) . But the 
map of IPhillips etai] ( 1 1998!) . with better positional accuracy and 
uv-coverage than in our measurement, reflects better the relative 
offset between the 6.7-GHz maser and the Hll region. 

The list of sources without a 9.9-GHz detection is given in 
Table |4l The first 7 columns are self-explanatory and include the 
source name, the date of observations, observed position and ve- 
locity range searched for maser emission as well as the synthesised 
beam size (given as the full width at half maximum) and its po- 
sition angle. The eighth column contains the noise level (la) in 
the dirty cubes searched for masers. This raw value does not take 
into account decorrelation due to weather, which we estimated us- 
ing the calibrator data (Table[T). As mentioned above, we adopted 
the largest estimate of the decorrelation factor obtained for the par- 
ticular observing slot. The cube noise level scaled with this factor 
is given in the ninth column of Table |4] as the estimate of the true 
la detection limit. As this limit is likely to be overestimated a bit 
for the majority of sources observed in a bad weather we give the 
mean temporal and angular separation from the appropriate sec- 
ondary calibrator in the last columns of Table |4] The smaller the 
separation, the lower the actual detection limit is expected to be. 
The uncertainties are given in the round brackets and expressed in 
the units of the least significant figure. 



3.1 Notes on selected targets 

In this section we discuss the sources where the 9.9-GHz masers 
were detected along with some non-detections. The latter are 
largely star-forming regions where an energetic interaction between 
outflow and molecular cloud has been reported in the literature, 
and therefore conditions favouring detection of the 9.9-GHz maser 
could be expected. 



3.1.1 G301. 14-0.23 (IRAS 12326-6245) non-detection 

This source is deeply embedded in a dense molecular cloud and 
harbours one of the most energetic and massive bipolar molec- 
ular outflows known to da te among objects of similar lumi- 
nosity tHenning et alj [2OO0I) . There is a weak (~1 Jy) class II 
methanol maser at 6 . 7 GHz with a p eak velocity of —39.8 km s~^ 
( ICaswelletal]|l995L ICaswellll2009l) . It is associated with one of 
the tw o compact Hll regions present in the source (Wal sh et al.l 
Il998h . Within the uncertainty of the position measurement this 
maser is co-located with the mainline OH maser, which has fea- 
tures spread over unusually large velocity range fro m —64 to 
—35 km s~ 1 dCaswell & Havneslll98llCaswelllll998l) . Voronkov 
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Table 3. The only four known masers at 9.9 GHz. In addition to 3 sources reported in this paper we have included G343. 12—0.06 jVoronkov et alj|2006l) 
for completeness. The uncertainties in the absolute positions given in the table are based on the formal uncertainty of the fit and the uncertainty of calibrator 
position. The commonly adopted systematic uncertainty (Icr) of the absolute positions measured with ATCA is C'S. All uncertainties in the table are given 
in the units of the least significant figure. The flux density uncertainty includes both the image noise (see Table|2) and the accuracy of the absolute flux scale 
calibration. The uncertainty of the peak velocity is the spectral resolution. 



Source 


Date of 






Determined position 


Peak 


Flux 


name 


observation 


1 


b 


^2000 


"52000 


velocity 


density 






(°) 


(°) 




(° ' ") 


(kms-i) 


(Jy) 


0331. 13-0.24 


16 Dec 2005 


331.131 


-0.245 


16:10:59.76 (5) 


-51:50:26.9 (4) 


-91.16 


1.92 (7) 




2 Jun 2008 






16:10:59.51 (3) 


-51:50:26.1 (3) 


-91.15 


3.1 (1) 




20 Jul 2008 






16:10:59.63 (3) 


-51:50:27.3 (3) 


-91.16 


2.6 (1) 


0343. 12-0.06 


16 Jun 2005 


343.121 


-0.064 


16:58:16.460(2) 


-42:52:25.73 (3) 


-31.56 


9.5 (3) 


W33-Met 


25 Mar 2006 


12.796 


-0.192 


18:14:10.897 (3) 


-17:55:58.42 (8) 


+32.72 


4.3 (1) 


O19.61-0.23 


25 Mar 2006 


19.608 


-0.232 


18:27:37.475 (2) 


-11:56:37.77 (8) 


+41.24 


3.3 (1) 



et al. (unpublished observations) have detected a class I methanol 
maser at 25 GHz (J=5 transition of the 25-GHz series) peaking at 
—36.0 km s~^ in this source. No 9.9-GHz maser emission has been 
detected in this source with the la flux density limit of 0.4 Jy. 



3.1.2 G305.21+0.21 possible detection 

The class II methanol maser at 6. 7 GHz in this sou rce was ob- 
served at high angular resolution bv' Norris et alj ( fl993h . who found 
two centres of activity separated by approxi mately 22", labelled 
G305A and G305B following the notation of I Walsh et all ilOOH) 
or 30 5.208+0.206 and 305.202+0.208 respectively according to 
ICasw ell (2009). There are two prominent Hll regions in the source 
located approximately 30" to the south-east of G305A and 15" to 
the west of G305B, so none of the maser sit es appear to be di- 
rectly associated w ith any continuu m emission i Phillips et aljl99^ : 
IWalsh et al.1 12007|) . According to 'Wals h et alj ilOOH) the region 
G305A is a hot core active in a number of molecular tracers, while 
the state of G305B, which is associated with bright 8 fim emission 
is not fully clear. 

This source was observed during the test 9.9-GHz survey 
without proper imaging. The spectrum contains a hint (at 1. 9a) 
of maser emission near —42.5 km s^^ with flux density around 
0.2 Jy. The regular su rvey was not sensitive enough to confirm 
the detection. However. I Voronkov et al.l | |2007|) reported a convinc- 
ing detection of the 104 GHz methanol maser in the source at the 
same velocity. The 104 and 9.9-GHz transitions belong to the same 
(J-l-l)_i — J_2 E transition series (J=10 and J=8, respectively) and 
are expected to show similar behaviour. Therefore, it is likely that 
the 9.9-GHz maser is present in this source below the detection 
threshold of the survey. The flux density of the 9.9-GHz emission is 
at least an order of magnitude lower than might have been expected 
on the basis of the only two other sources k nown to have both 104 
and 9.9-GHz masers. Both G343. 12 -0.06 j Voronkov et"^l2006l) 
and W33-Met ( I Voronkov et al.ll2007h have flux densities at 104- 
GHz comparable with the 9.9-GHz ones reported in this paper (see 
Table[3}. This most likely implies that markedly different physical 
conditions exist in this particular source. One may s peculate that 
this m ight be related to another peculiarity found by IWalsh et al.l 
the high-mass star-forming region associated with G305B 
is bright in infrared, old according to the chemistry, but has no Hll 
region sufficiently developed to be detectable. The expansion of 
an Hll region could be hindered because of the dense environment 
and/ or by continuing infall o f the matter. 

IVoronkov et al.l ( l2005al) detected a bright class I methanol 



maser at 25 GHz (J=5 transition of the 25-GHz series) at 
—42.4 km s~^, although without imaging and accurate flux den- 
sity calibration. These masers have been shown to originate at the 
same spatial location as the 9.9-GHz maser in at least one source 
jVoronkov et alj2 006V An accurate position for the 25-GHz maser 
(J=4 transition with the peak velocity of —42.1 km s^^ observed 
with the velocity r esolution of 0.75 km s~^) was determined by 
IWalsh et"ail ( l2007h . but the reported flux density of the maser is 
much lower, due to insufficient spectral resulution. The 25-GHz 
maser was found to be 3" to the east from the position of the 
6.7 GH z maser 305.202+ 0.208 (G305B region). It is worth not- 
ing that lWalsh et all l l2007h also detected 25-GHz emission (several 
transitions of the series) towards two other locations in this source: 
near G305A and the south-eastern Hll region. However, the radial 
velocity and a clear maser appearance of the feature near G305B 
allows an unambiguous identification of the suspected 9.9-GHz 
maser with this position, near the 6.7-GHz maser 305.202+0.208. 



3.1.3 G3 3 1.1 3 -0.24 detection 



The 6 .7-GHz masers in this source were mapped bv lPhillips et al.l 

and found to be located near the edge of a 7x 10 arcsec 
Hll region (see also Fig.|2^). According to .Goedhart et al.l ( l2004h 
this 6.7-GHz maser shows periodic flares with period around 500 
days. The OH maser in the source is coincident with the cluster of 
the 6.7-GHz maser sp ots within the uncertainty of the a bsolute po- 
sition measurement (Philli ps et al]| 19981 : ICaswell|[l998 ). The 9.9- 
GHz maser is located close to the cluster of class II masers, but at 
a clearly distinct position. It is also seen projected onto the Hll re- 
gion but located near the edge. There is a protrusion in the image 
of the Hll region near the location of the 9.9-GHz maser (Fig.|2^) 
indicating a possible interaction between expanding Hll region and 
the surrounding material. The brightest 4.5-p.m source within the 
boundaries of this HlI region is located between the 6.7 and 9.9- 
GHz masers (Fig. |2^). Two weaker 4.5-/im sources surrounded 
by an extended emission are located near the northern tip of the 
Hll region (as shown by an ellipse in Fig. |2^). This latter object 
has an excess of the 4.5- fim emission and is listed as an EGO by 
ICyganowski et al.l ( 1200 8l) despite being only a few seconds of arc 
across. The exact physical interpretation of this source is uncertain 
at present. However, considering its location it seems unlikely that 
this EGO could trace part of the shock front responsible for the 

9.9- GHz maser. 

iDe Buizer et aD ( l2009l) mapped the SiO (2-1) emission to- 
wards G331.13— 0.24 and detected an elongated SiO source with 
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Table 4. Sources without a 9.9-GHz detection. Flux density limit is Icr. G329. 031-0. 19 and G329. 029-0. 20 represent two close positions within the same 
source. 



Source 


Date of 


Observed position 


Velocity 


Synthesised beam 


Icr image 


lo" flux 


Separation 


from calibrator 


name 


observations 


C12000 


^^2000 


range 


FWHM 


pa 


noise level 


limit 


angular 


temporal 








^/i m 


(° ' ") 


(km s-i) 


(arcsec) 


(deg) 


(mjy beam~^) 


(mJy) 


(deg) 


(min) 




G301. 14-0.23 


15 Dec 2005 


12:35:34 


-63:02:38 


-97, -1 


1.1x0.8 


39 


76 


430 


10.5 


6.6 


(1) 


G305.21-I-0.21 


15 Dec 2005 


13:11:09 


-62:34:41 


-97, -1 


1.1x0.8 


33 


82 


470 


14.5 


6.9 


(2) 


G305.25-I-0.25 


15 Dec 2005 


13:11:33 


-62:32:03 


-97, -1 


1.3x0.7 


34 


76 


430 


14.6 


6.8 


(2) 


G305.36-I-0.20 


15 Dec 2005 


13:12:34 


-62:33:25 


-97, -1 


1.3x0.7 


35 


81 


460 


14.7 


17.2 


(1) 


G326.475-I-0.70 


15 Dec 2005 


15:43:17 


-54:07:13 


-94, +2 


1.3x0.7 


24 


76 


430 


6.6 


6.5 


(1) 


G326.641+0.61 


15 Dec 2005 


15:44:33 


-54:05:29 


-94, +2 


1.3x0.8 


26 


69 


390 


6.5 


6.4 


(1) 


G327.392+0.19 


16 Dec 2005 


15:50:18 


-53:57:06 


-124, -28 


1.6x0.9 


20 


68 


460 


6.1 


6 


(1) 


G326.859-0.67 


16 Dec 2005 


15:51:14 


-54:58:05 


-125, -29 


1.6x1.0 


9 


89 


610 


5.2 


5.4 


(2) 


G327.618-0.il 


16 Dec 2005 


15:52:50 


-54:03:01 


-124, -28 


1.5x1.0 


14 


64 


440 


5.8 


5.2 


(1) 


G327.29-0.58 


15 Dec 2005 


15:53:06 


-54:37:06 


-95, -fl 


1.3x0.7 


25 


62 


350 


5.3 


6 


(1) 


G328.81-I-0.64 


15 Dec 2005 


15:55:48 


-52:42:48 


-94, +2 


1.3x0.8 


32 


59 


340 


6.8 


16 


(3) 


G328.237-0.54 


17 Dec 2005 


15:57:58 


-53:59:23 


-94, +2 


1.2x0.8 


-16 


52 


98 


5.5 


5.0 


(1) 


G329.469-I-0.50 


16 Dec 2005 


15:59:41 


-52:23:28 


-124, -28 


1.5x1.0 


14 


61 


410 


6.9 


5.3 


(1) 


G329.031-0.19 


17 Dec 2005 


16:00:30 


-53:12:27 


-94, +2 


1.2x0.8 


-14 


52 


98 


6.1 


4.9 


(1) 


G329.029-0.20 


17 Dec 2005 


16:00:32 


-53:12:50 


-94, +2 


1.2x0.9 


-12 


52 


98 


6.1 


4.9 


(1) 


G329.066-0.30 


17 Dec 2005 


16:01:10 


-53:16:03 


-94, +2 


1.2x0.8 


-11 


53 


100 


6.0 


4.9 


(1) 


G329.183-0.31 


17 Dec 2005 


16:01:47 


-53:11:44 


-94, +2 


1.1x0.9 


-3 


52 


98 


6.0 


5 


(1) 


G33 1.442-0. 18 


16 Dec 2005 


16:12:12 


-51:35:10 


-124, -28 


1.8x0.9 


16 


51 


350 


7.3 


7 


(3) 


G33 1.342-0.34 


16 Dec 2005 


16:12:26 


-51:46:17 


-124, -28 


1.5x0.9 


36 


61 


410 


7.1 


5.3 


(1) 


G332.295-0.09 


17 Dec 2005 


16:15:45 


-50:55:54 


-94, +2 


1.1x0.9 


34 


52 


98 


7.9 


5.3 


(8) 


G332.604-0.16 


17 Dec 2005 


16:17:29 


-50:46:13 


-93, -1-3 


1.2x0.8 


-7 


53 


100 


8.0 


5.0 


(1) 


G333.029-0.06 


17 Dec 2005 


16:18:57 


-50:23:54 


-93, -f 3 


1.2x0.9 


-9 


52 


98 


8.4 


5.0 


(1) 


G333.315-I-0.10 


17 Dec 2005 


16:19:29 


-50:04:41 


-93, 4-3 


1.2x0.9 


-8 


52 


98 


8.7 


5.0 


(1) 


G333. 184-0.09 


16 Dec 2005 


16:19:46 


-50:18:35 


-124, -28 


1.7x1.0 


29 


57 


390 


8.5 


5.0 


(5) 


G333. 163-0. 10 


25 Mar 2006 


16:19:43 


-50:19:53 


-136, -40 


1.3x0.7 


5 


61 


76 


4.9 


5.3 


(1) 


G333.23-0.05 


16 Dec 2005 


16:19:48 


-50:15:02 


-124, -28 


1.8x1.0 


20 


54 


370 


8.6 


5 


(1) 


G333.121-0.43 


17 Dec 2005 


16:21:00 


-50:35:52 


-94, -1-3 


1.1x0.9 


34 


49 


92 


8.2 


5.0 


(1) 


G333.562-0.02 


17 Dec 2005 


16:21:09 


-49:59:48 


-93, -1-3 


1.2x0.9 


-3 


52 


98 


8.8 


5 


(1) 


G332.942-0.69 


17 Dec 2005 


16:21:19 


-50:54:10 


-94, +2 


1.1x0.9 


29 


49 


92 


7.9 


5.3 


(5) 


G333.466-0.16 


25 Mar 2006 


16:21:20 


-50:09:49 


-105, -9 


1.3x0.7 


8 


60 


75 


4.6 


5.2 


(1) 


G332.963-0.68 


25 Mar 2006 


16:21:23 


-50:52:59 


-106, -10 


1.3x0.7 


9 


61 


76 


4.6 


13.9 


(1) 


G333.130-0.56 


25 Mar 2006 


16:21:36 


-50:40:51 


-106, -10 


1.4x0.7 


12 


61 


76 


4.5 


5.2 


(1) 


G335.060-0.42 


26 Mar 2006 


16:29:23 


-49:12:27 


-75, -1-22 


1.5x0.7 


19 


57 


71 


3.7 


5.4 


(1) 


G337.40-0.41 


25 Mar 2006 


16:38:50 


-47:28:18 


-103, -7 


1.4x0.7 


9 


61 


76 


3.8 


5.1 


(1) 


G338.92+0.56 


16 Dec 2005 


16:40:32 


-45:41:53 


-123, -27 


2.0x0.9 


8 


50 


340 


13.6 


5.7 


(3) 


G338.92-0.06 


25 Mar 2006 


16:43:16 


-46:05:42 


-102, -6 


1.4x0.7 


9 


61 


76 


4.8 


13.8 


(1) 


G339.62-0.12 


26 Mar 2006 


16:46:05 


-45:36:43 


-72, -1-24 


1.6x0.7 


17 


57 


71 


5.2 


5.6 


(1) 


G345.01-I-1.79 


26 Mar 2006 


16:56:49 


-40:14:20 


-69, -1-27 


1.7x0.7 


16 


58 


73 


10.6 


5.9 


(3) 


G345.00-0.22 


26 Mar 2006 


17:05:10 


-41:29:04 


-69, 4-27 


1.7x0.7 


16 


58 


73 


9.6 


15.2 


(3) 


G351. 16-0.70 


26 Mar 2006 


17:19:57 


-35:57:46 


-66, 4-30 


1.9x0.7 


15 


58 


73 


9.1 


6.0 


(6) 


G351.24-I-0.67 


25 Mar 2006 


17:20:16 


-35:54:58 


-24, -f46 


1.4x1.0 


30 


62 


78 


9.0 


5.5 


(3) 


G35 1.4 1-1-0.64 


26 Mar 2006 


17:20:53 


-35:47:02 


-66, -1-30 


1.9x0.7 


15 


58 


73 


8.9 


5.9 


(4) 


G351.78-0.54 


26 Mar 2006 


17:26:43 


-36:09:18 


-66, -1-30 


1.8x0.7 


16 


57 


71 


9.6 


5 


(1) 


GO.54-0.85 


25 Mar 2006 


17:50:15 


-28:54:51 


-32, -1-64 


1.7x1.0 


23 


60 


75 


8.3 


14.5 


(5) 


G9.621-I-0.196 


26 Mar 2006 


18:06:15 


-20:31:32 


-60, -1-37 


3.0x0.7 


9 


58 


73 


13.7 


16 


(1) 


G18.34-I-1.77 


25 Mar 2006 


18:17:57 


-12:07:22 


-28, -1-68 


3.7x1.1 


7 


61 


76 


11.0 


13 


(3) 



velocity gradient centred at the position of the 6.7-GHz masers 
with the axis roughly parallel with the E ast- West distribution of the 
6.7-GHz maser spots. IPe Buizeil (l2003h reported the detection of a 
single knot of H2 2.12-/im emission. This is located at the eastern 
edge of the red-shifted SiO lobe (approximately 8" from the 6.7- 
GHz maser). Both these facts reveal pronounced shock activity in 
the region. However, the present data do not allow us to rule out the 
association of the 9.9-GHz maser with any certain structure seen 
in shock-tracing molecular emission as the spatial resolution of the 
SiO data is rather coarse and comparable to the separation between 
the 9.9- and 6.7-GHz masers. 



3.1.4 G338.92-0.06. G339.62-0.12, G9.62+0.20 
non-detection 

These sources were in cluded in the target li st as perio dically vari- 
able 6 .7-GHz masers jGoedhart etar]|2004 . although [sivsh et al.l 
l l 19941) detected no class I maser emission at 44 GHz towards ei- 
ther of them. iKurtz et 5] ( l2004h found weak 44-GHz maser emis- 
si on towards G9.62-I- 0.20, which was below the detection threshold 
of lSlvsh etal] ( ll994l) . Despite the fact that we observed all three 
sources in a good weather providing a 1-cr flux density limit of 
around 70-80 mJy (see Table |4ll, no emission at 9.9 GHz was de- 
tected. 
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3.1.5 G343. 12-0.06 (IRAS 16547-4247) detection 



G331, 13-0,24 




G343. 12-0.06 



-a 

X 




Vlsr^ km s" 

W33-Met 



-i-> 

m 
QJ 

X 



1 1 1 1 1 r 




36 



Vlsr. km 
G19.61-0.23 




Vlsr. km s"' 



Figure 1. Spectra of the masers detected at 9.9 GHz. The uncertainty in 
radial velocity due to the rest frequency uncertainty is less than the spec- 
tral r esolution. Fo r completeness, we show the spectrum of G343. 12—0.06 
ffrom lVoronkov et al. 2006). Note, the spectrum of this source was obtained 
with four times better spectral resolution, which is comparable with the rest 
frequency uncertainty. 



The 9.9-GHz maser was discovered in this source during the pi- 
lot survey. It has bee n investigated in detail in a separate paper 
dVoronkov etai]|2006ll and is included here for completeness. The 
maser is associated with the brightest knot of th e H2 emission 
terooks et al.l2003h tracing an outflow (see Fig.[3h: IVoronkov et al I 
llOOd) . This particular location is the only place where the 25 
and 104-GHz masers were detected in this source, in contrast to 
the widespre ad class I methanol m asers detected in a n umber of 
maser spots jVoronkov et"ai] l2006h . [Garay et al] ( |2003|) first de- 
tected cm-wavelength continuum emission in the source compris- 
ing a brig hter central object and two weaker satellites (further re- 
solved bv lRodrfguez et al]|2005t) . The central object and the satel- 
lites were interpreted as thermal (free-free) emission from a ra- 
dio jet and internal working surfaces of the jet, r espectively. The 
12-mm image of these continuum sources from IVoronkov et al.l 
( I2OO6I) is shown in Fig.[3}i along with the 2.\2-^m molecular hy- 
drog en image jS rooks et al. 2003) and the 9.9-GHz maser posi- 
tion. Forster & Caswelll ([1989) detected mainline OH masers at a 
position close to the central continuum source (Fig. [3ji). No 6.7- 
GHz emission has been detected with an upper limit of 0.3 Jy 
dCaswell et aljl995l : lwaish et aljl998l) . 



3.1.6 G351. 16-0.70, G351.24+0.67, G351.41+0.64 
non-detection 

These are the positions within NGC 6334 complex, one of the 
most prominent sites of massive-star formation. The y correspond 
to so urces V, IV and I/I(N), respectively (see e.g. iMunos et al.l 
I2OO7I and references therein). The massive twin cores NGC 
63341 and I(N) have been studied for more than two decades 
and are believed to be at different evolutionary stages (see 
ISeuther et al.ll2008l , and references therein). These cores are well 
known to host unusua l methanol masers of both classes (e.g ., 
lMenten&Batri3ll989l : lEUingsen et al] |2004 ICragg etalj 12004) . 
There are class II methan ol masers at 6.7 GHz associated with both 
cores dWalsh et alj[T99^ . but the maser in the core I is one of the 
brightest known 6.7-GHz masers of more than 3000 Jy. In contrast, 
the n orthern core, I(N), i s remarkable in its class I maser emission 
(e.g.. iBrogan et alj|2009h . Although both sites emit at 25 GHz, it 
is likely to be a maser in the I(N), while source I is responsible 
for quasi-thermal emission jMenten & Batrlall989l : lHaschick et"ail 
119901 : iBeuther et ar]|2005h . No 9.9-GHz maser emission has been 
detected towards either of these positions with the Icr flux density 
limit of 70-80mJy. 



3.1.7 W33-Met (G12.80-0.19) detection 

The W33 region, which hosts a cluster of young massive stars, is 
renowned for its strong radio continuum emission and a complex 
molecular environment jGoss et al 1 19781 : iGoldsmith&Maoll 19831 : 
Ho et al. 1986). The fir st detection of a meth anol maser towards this 
source was reported bv lMenten et al. I|l986). They observed a num- 
ber of class I maser transitions of the 25-GHz series and revealed a 
positional offset of the masers with respect to the dominant contin- 
uum source. The 44-GHz methanol masers were later found in the 
vicinity of this offset position by Haschick et al. (1990), who in- 
troduced the notation W33-Met to distinguish methanol maser site 
from the main continuum source in the region. Before we started 
the pr oject, this was the only source known to have a 9.9-GHz 
maser Jsivsh etalJ[T993l) . Although our measurement gives the 
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Figure 2. Radio continuum emission at 8.6 GHz (contours) in the vicinity of masers detected at 9.9 GHz overlaid on t op of tlie Spitzer IRAC 4.5-/xm image 
(grayscale). The environment of G343. 12—0.06 is shown in Fig. [5^ and discussed in a separate paper I Voronkov et al . 2006). Position of the 6.7-GHz and 
9.9-GH z masers is shown by the open squares and crosses, respectively, (a) G331. 13—0.24. The continuum image was obtained using the data o f lPhillips et ai] 
<1998l) . to preserve high accuracy of the relative position of the 6.7-GHz maser with respect to the Hll region. Contours are 1, 2, 4, 8, 16, 32, 64 and 90% of 
43 mj y beam~ ^ . The crosses r epresent 3 epochs of the measurement at 9.9 GHz, rather than individual maser spots. The ellipse represents the location of an 
EGO fcvganowski et alj2008l) . (b) W33-Met (G12.80— 0. 19). The 4.5-/im flux density has been multiplied by 10 inside the square in the bottom right corner 
of the image. Contours are 5, 10, 20, 30, 40, 50, 60, 70, 80 and 90% of 1.8 Jy beam^i. (c) G19.61-0.23. Contours are 5, 10, 20, 30, 40, 50, 60, 70, 80 and 
90% of 0.6 Jy beam-l. 



peak flux density as 4.3 Jy, ISlvshetal.l l [l993h reported just 0.8 Jy. 
Partia lly, this is the resul t of the 2.5 times lower spectral resolution 
m the lSlvshetaP i ll 9931) data, as the line is quite narrow (Fig. [T]l. 
However, the spectral resolution alone cannot account for the more 
than factor of 4 difference in the flux density. It is no t clear how 
accurate the absolute flux scale calibration was for the lSlvsh et al.l 
( 1 19931) single dish data, especially given the strong cm-wavelength 
continuum in the source. Therefore, we do not draw any conclu- 
sions from t his discrepancy, alt hough temporal variability cannot 
be excluded. ISlvsh et al.l ( [l993h reported the peak velocity of the 
maser at 33.5 km s~^, which is notably different from 32.72 km s~^ 
obtained in our measurement (Table [3]( even taking into consider- 
ation the spectral resolution achieved in both experiments. How- 
ever, this discrepancy can b e explained by the less accurate value 
of the rest frequency used bv lSlvsh et al.l i 19931) to calculate the ra- 
dial velocities. The difference in adopted rest frequencies accounts 
for 0.85 km s~^ of the velocity offset, bringing the radial velocity 
into agreement within the spectral resolution. The radial velocity 
of the 9.9-GHz maser is als o in agreement with the velocity of the 
25-GHz masers report ed by Menten et alj ( 19861) . provided the rest 
frequency adopted by iMenten et al.l jl98( ) is more accurate than 
the spectral resolution of their measurement. 

The continuum source has a complex structure reflecting the 
presence of multiple young stars (Fig|2j3). The mo rphology is in 
general agreement with the earlier results of IHo et alj l il986l) , but 
the 1.5-GHz emission (20cm) is extended more to the west than the 
emission at higher frequencies. This is most obvious if one com- 
pares the association between the (1,1) inversion transiti on of am- 
mon ia and the rim of the 1.5-GHz emission depicted by K eto & IJol 
( Il989 ) in their Fig. 2 with Fig.jSji in this paper, where the same am- 
monia distribution from , Keto & Ho tl98ft) is overlaid on top of the 
8.6-GHz continuum from our observations. Fig.[3}3 clearly shows a 
gap of around 15" between the 8.6-GHz continuum source and the 
ridge of ammonia emission. As concluded by H o et a l. ( 1986), the 
smooth emission filling the gap is most likely resolved at higher 
frequencies. According to the interpretation of iKeto & Hoi jl989l) 
the ridge of ammonia emission is due to shocks caused by interac- 
tion of the ionised gas with the neutral component. The 9.9-GHz 



maser is located near the western clump of the ammonia emission 
(Fig. 1315). This position is in agreement with t he position of the 
25-GHz masers reported bv lMenten et al.l ( Il986l) within the uncer- 
tainty of their measurement. 

The Spitzer IRAC image revealed a number of 4.5-^m sources 
associated with some of the continuum peaks (see Fig|2j5). There is 
an extended 4.5-/im emission located largely near the western edge 
of the 8.6-GHz continuum source in the area where the continuum 
emission has been detected at 1.5 GHz only. The inspection of other 
IRAC bands revealed no significant excess of the 4.5-/im emission, 
w hich is in agreement with this object not being listed as an EGO 
bv lCvganowski et al] ( |2008|) . However, it is worth noting that there 
is a faint compact 4.5-/im source located just 3" to the west of the 
9.9-GHz maser position. This infrared source is detected in the 4.5- 
^m IRAC ban d only. 

Although iMentenI ( Il99ll) reported the detection of the 6.7- 
GHz maser in W33-Met, a compari son of the ve locity range with 
the source G12.91— 0.26 reported bv lCaswell et a l. ( 1995) suggests 
possible confusion. We processed archival ATCA data (project 
code C415) using the standard reduction technique and found two 
6.7-GHz masers G12.68-0.18 (Q2000 = 18''13'"54';74±0;01, 
^2000 = -18°01'46'.'l±0'.'l) and G12.91-0.26 (Q2000 = 
18'n4™39:52±0;02, <52ooo = -17°51'59'.'l±0'.'2), situated 7!0 
and 7^9 away from the 9.9-GHz methanol maser position, respec- 
tively. The la rms n oise level at the position of the 9.9-GHz maser 
at 6.7 GHz is 0. 1 Jy. ICaswelll ( |2009|) provided an independent po- 
sition measurement of the confusing 6.7-GHz masers which agrees 
well with the values given above. 



3.1.8 G18.34+1.77 (IRAS 18151-1208) non-detection 

Although the widespread class I maser transitions (e.g. at 44 GHz) 
are not atypical in this source, it has one of the brightest 25- 
GHz methanol masers found outside Orion and NGC 6334 (the 
J=5 maser peaks at +29.7 km s^^ and has a flux den sity exceed- 
ing 15 Jy, Voronkov et al., unpublished observations). iDavis et al.l 
(2004) detected two H2 outflows in this source which could be re- 
lated to the 25-GHz masers. High angular resolution observations 
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are required to confirm this. No 9.9-GHz maser emission has been 
detected in this source with the la flux density hmit of 76 mjy, 
which imphes that the flux density of the 25-GHz maser exceeds 
that of the 9.9-GHz maser by at least two orders of magnitude. 



3.1.9 G19.61 -0.23 detection 

This region hosts a cluster of young high-mass stars and is some- 
what similar in appearance to W33-Met at cm-wavelength (Fig.[2j3 
and c) , albeit with a sm aller angular separation between compo- 
nents. iFuruva et al.l ( 12005 ) identified nine ultra-compact Hll re- 
gions and a hot core squeezed between two of them designated A 
and C. The 8.6-GHz continuum image (Fig. [2}:) obtained in ou r 
survey qualitatively agrees with the results o f iFuruva et all ilOO^ . 
although the elongated beam does not allow us to disentangle close 
components easily. A number of continuum peaks coincide with 
the 4.5-/im emission from the Spitzer IRAC image (see Fig|2};). 
The brightest 4.5-^m sour ce is associated with the most prominent 
Hll region (labelled A in Furuva et al. 2005). This source is not 
listed as an EGO by .Cyganowski et al. (2008). 

The hot core is responsible for emission in a number of molec- 
ular species dFuruva et al]|2005l ; IWu et alj|2009l) and is associated 
with the brightest ammonia clump (nam ed middle clump) found by 
iGarav et al.rjl998l) . lGarav et alj ( 1199 8') identified two more ammo- 
nia clumps: the northern clump seen in absorption and the elon- 
gated south-western clump seen in emission. The rela tive position 
of th e 8.6-GHz continuum, ammonia emission (from lOarav et al.l 
1 1998 absorption source is not shown) and the 9.9-GHz maser is 
shown in Fig.[3j;. The 9.9-GHz maser is l ocated near the elongated 
south-western clump, which iGarav et alj j 19981) interpreted as aris- 
ing from the compressed and heated gas behind the shock front 
driven by the expansion of the largest Hll region in the source. 
I Kurtz et a l. ( 2004) observed this source at 44 GHz and revealed two 
maser spots. The brightest 44-GHz spot is co-located both in posi- 
tion and velocity with the 9.9-GHz maser within the combined error 
of the measurements. The weaker feature peaks at —46.5 km s~^. 
It corresponds to a maser spot located 7'.' 9 to the North-East from 
the brightest spot (and the 9.9-GHz maser spot) roughly along the 
am monia filament. 

ICaswell et al] d 19951) detected a weak 6.7-GHz methanol 
maser with a peak flux density of 0.4 Jy in this source. This maser 
is located 7'.'6 away from the position of the 9.9-GHz maser, at 
02000 = 18''27'"37^993±ff:006, S2000 = -11°56'37'.'6±0'.'4 
(Caswell, unpublished observations). The position suggests an as- 
sociation of this maser wit h the hot core and the middle am- 
monia clump (see Fig. [3};). iGarav et"ail jl985l) found a number 
of OH maser spots clustered near the 6.7-GHz maser position 
and also associated wi th the middle ammonia clump (see also 
iForster & Caswellll 19891) . 



variable masers known to date dGoedhart et"ai]|2004h . which has 
the longest period in the sample of around 500 days. The coexis- 
tence of 6.7 and 9.9-GHz masers in this source is a unique situation 
with great potential to shed light on the nature of periodic flares, be- 
cause these two transitions are predicted to react in opposite ways 
to the changes in pumping. The pumping mechanisms of these two 
maser transitions are in conflict as they belong to different classes. 
Hence, if the flares occur due to a pumping variation of any nature, 
the intensity variations in the two transitions are expected to anti- 
correlate. From the other hand, both masers are seen projected to- 
wards the same Hll region. This means that this Hll region provides 
the seed radiation which both masers amplify. If the variations are 
due to changes in the seed radiation, intensity variations in the two 
transitions will correlate. Thus simple monitoring of the 9.9-GHz 
maser may unravel a vital aspect of the flare mechanism. We have 
started such a project and will report the results elsewhere. 

The main complication arises from the fact that class I and 
class II masers are not co-located. In this and some other peri- 
odically variable sources, flares of the different 6.7-GHz spectral 
features, which are presumably also separated in space, are de- 
layed with respect to each other, sometimes by up to a few weeks 
(Goedhart et al. 2004). Although the velocity of the 9.9-GHz maser 
is close to the velocity of one of the groups of the 6.7 GHz maser 
lines, the angular separation is an order of magnitude greater than 
the size of the 6.7-GHz cluster. Therefore, a longer delay of the or- 
der of months is quite possible. The 9.9-GHz flux measurements 
presented in this paper all correspond to the descending part of the 
6.7-GHz light curve, approximately 20, 4 and 10 weeks after the 
maximum (Goedhart et al., unpublished recent monitoring data), 
respectively for each epoch listed in TablejS] Although this is con- 
sistent with the steady fall of the 9.9-GHz flux density along with 
the 6.7-GHz flux density, a longer time series is required to obtain 
a definitive answer. The long period of around 5 00 days as well 
as the non-sinusoidal shape of the light curve (see ^ Goedhart et akl 
2004) will help to distinguish between a possible delay and anti- 
correlation in this particular source. 

Monitoring of the masers has a number of advantages over 
monitoring the continuum emission directly. Firstly, the masers am- 
plify the continuum and so essentially undetectable variations of 
the continuum flux can cause significant variations in the maser 
intensity. Secondly, the continuum source has a complex shape in 
G33I.I3— 0.24 making it very difficult to undertake accurate flux 
density monitoring as it is hard to reproduce exactly the same uv- 
coverage in all epochs. In contrast, each spectral feature of the 
methanol maser corresponds to a point source for the baselines of 
the order of a few kilometres long, making flux monitoring a much 
easier task. It is worth noting that there is an agreement to better 
than 2% (which is likely to be the limit of the absolute flux density 
calibration) between the flux density of the Hll region measured by 
IPhillips et al., (.1998.) and that obtained in our study. 



4 DISCUSSION 

4.1 Temporal variability of G331.13 0.24 

Inspection of the 3 epochs listed in Table [3] suggests possible tem- 
poral variability of the 9.9-GHz maser in G331.13— 0.24, with the 
caveat that the only measurement which stands out significantly 
was taken under imperfect weather conditions. The variability of 
the class I methanol masers is not well understood at the moment. 
However, the class II maser at 6.7 GHz in this source is known 
to be somewhat exceptional. This is one of the seven periodically 



4.2 Coexistence with class II methanol masers at 6.7 GHz 

As was mentioned above, our target list comprised largely known 
class I methanol masers. The majority of these class I masers were 
found towards known 6.7-GHz masers (class II), as they are often 
located within the beam of the single dish used for class I mea- 
surements (e.g. Slvsh et al. 1994; Ellingsen 2005). The remaining 
sources were either well known star-forming regions (e.g. posi- 
tions within the N GC6 334 compl ex) or were found towards known 
H2O masers (e. g. iBach iller et al. ) ll990l) . Hll regions or main-line 
OH masers (e.g. lSlvsh et al.lI994h . Despite being selected by other 
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Figure 3. Association of the 9.9-GHz masers witli outflow and hot molecular gas in the vicinity of Hll regions. Position of the 9.9-GHz masers is shown by 
crosses and the 6.7- GHz maser in just G19.61 — 0.23 by an open squar e, (a) G343. 12—0.06. The maser position is overlaid on top of the H2 2.12-/im emission 
ifirooks et al.ll20O3h and the 12mm continuum iVoronkov et al.l2006ll showing an a ssociation with ou tflow. Contours are 15, 50 and 90% of 13 mjy beam~^ . 

(b) W33-Met (G12.80— 0. 19). The grayscale shows the (1,1) ammonia emiss ion jKeto & Holl 19891) . Contours show the 8.6-GHz continuum as in Fig|2j>, 

(c) G19.61— 0.23. The grayscale represents the (2,2) ammonia emission from iGaravetalTj 19981) . The northern molecular clump seen in absorption is not 
shown in this figure. Contours show the 8.6-GHz continuum as in Fig|2];, 



criteria many of the sources from the second group also have a 
6.7-GHz methanol maser in the vicinity. Therefore, only one tar- 
get in our sample, W33-Met, has a class I methanol maser with- 
out a 6.7-GHz counter part. Another such source is G343.12— 0.06 
jVoronkov et"ai]|2006h . 

In total, out of four 9.9-GHz methanol masers currently 
known, two do not have a 6.7-GHz counterpart detected, and one 
source, G19.61— 0.23, has a very weak 6.7-GHz maser. The fourth 
source, G331.13— 0.24, has both masers projected onto an HlI re- 
gion, a situation favouring detection of low gain masers. These re- 
sults suggest that the 9.9-GHz masers tend to be associated with 
sources with class I masers without strong class II counterparts. 
A proper statistical test is not yet possible due to a small num- 
ber of known 9.9-GHz masers. However, a number of previous 



studies have already investigated relations between the widespread 
class I mase rs (e.g. at 44 or 95 GHz) and the class II masers a t 
6.7 GHz (e.g. lSlvsh et al.ll999l : lEllingsenl2005h . lSlvsh et alj ( ll999h 
were motivated by the conflicting pumping mechanisms of these 
masers, although the present data suggest that the class I and class II 
mas ers are rarely seen at the same location at hi gh spatial resolution 
(e.g. lKurtz et alj2004l : lcv"ganowski et al.l2009l) . The exceptions are 
probably due to the projection effect, although unusual pumping 
under common conditions (e.g. Voronkov et al. 2005b) cannot be 
excluded. Therefore, any tendency (e.g. for bright masers of one 
class to avoid the bright masers of another class), if present, is likely 
to be related to the general trends in the distribution of the physical 
parameters affecting the pumping across the star-forming region, 
rather than the conditions at any single location in the source. These 



New 9.9-GHz methanol masers 11 



trends can be idiosyncratic to some extent due to the complex struc- 
ture of high-mass star-forming regions (e.g. one region could be 
more embedded in the parent molecular cloud than another). How- 
ever, they can also change during the course of evolution, for exam- 
ple when an outflow turns on or an Hll region starts to expand. It is 
worth not i ng tha t using a statistically complete sample of sources, 
lEningseij ( l2003) found no evidence for the original anti-correlation 
between flux densities reported bv lSlvsh et al.l lll999'). However, the 
9.9-GHz masers require a more restrictive range of physical condi- 
tions (e.g temperature and density) to appear than the wide-spread 
class I methanol masers at 44 and 95 GHz. It is therefore possi- 
ble that the anti-correlation with the 6.7-GHz masers is more pro- 
nounced for these 9.9-GHz masers than for the widespread 44 and 
95-GHz ones. 



4.3 Association with expanding Hii region 

Out of 48 sources searched for the 9.9 GHz emission in this study 
(note that G329.031-0.19 and G329.029-0.20 are counted as a sin- 
gle source), 35 sources have 8.6-GHz continuum sources, presum- 
ably Hll regions, within the primary beam (half power width is 
about 5 arcmin) and 13 sources do not (Icr rms noise level not 
corrected for decorrelation and primary beam attenuation is about 
0.3 mjy). Three of these continuum sources correspond to periodi- 
cally variable 6.7-GHz masers and two of them do not have class I 
methanol masers detected so far in any transition. Even allowing 
for a few confused cases due to a large beam (most known class I 
masers have single-dish positions), our target sample, comprising 
largely known class I masers at 44 and 95 GHz, appears to be bi- 
ased to some extent towards the star-formation sites with Hll re- 
gions. The present data do not yet allow us to establish reliably 
whether the common occurrence for class I masers to have an Hll 
region in the vicinity has physical significance or is just a selection 
effect. In contrast, it is well established that only a small fraction of 
class II methanol masers at 6.7 GH z are associated with detectable 
Hll regions (e.g.. lWalsh etal]il998h . 

All four currently known 9.9-GHz maser sources have a cm- 
wavelength continuum detected in the vicinity of the maser spot. 
Three masers found in the regular survey are located in close prox- 
imity to well-developed Hll regions (Fig. |2j. While the remaining 
source, G343.I2— 0.06, is reno wned for thermal emission from the 
jet dr iving an outflow (Fig.[3b; lBrooks et aI.ll2003l : IVoronkov et al.l 
l2006h . As mentioned in section [3TT] in some cases the molecular 
emission could arise due to interaction between ionised and neutral 
material. Such an interpretation has been suggested on the basis of 
kinematics an d morphology for at least two source s, W33-Met and 
GI9.6I-0.23 jKeto & Holl989l : lGarav et al.ll998l) . The expansion 
of an Hll region drives shocks into the surrounding molecular cloud 
heating and compressing the gas. It is this change in the physi- 
cal conditions which is most likely traced by the ammonia emis- 
sion (Fig.[3l3c). In both W33-Met and G19.61-0.23, the 9.9-GHz 
maser is located just beyond the edge of the ammonia source, fur- 
ther away from the Hll region which is believed to be undergoing 
expansion. The separation between the maser and the nearest am- 
monia emission is of the order of a few seconds of arc. Although 
the magnitude of this separation could probably be explained by 
the combined accuracy of the absolute positions, the picture agrees 
well with the basic expectation that a more fragile molecule like 
methanol should be located predominantly further away from the 
exci ting source than ammon ia. 

iKirsanova et 1a\ j2009l) modelled the expansion of an Hll re- 
gion surrounded by molecular gas and demonstrated formation of 



a transition layer characterised by high abundance of molecules in 
the gas phase. The location and thickness of this transition layer 
vary from molecule to molecule and also depend on the initial tem- 
perature and density of the qui escent gas. Although th e methanol 
chemistry was not modelled by IKirsanova et all ( 120091) . the physi- 
cal conditions in the transition layer, such as temperature and den- 
sity, simulated in these models could cle arly be in the range su itable 
for the formation of class I masers (e.g..'Voron kov et ahl 200 5b). 

The geometry of G331.13— 0.24 and the data on the thermal 
molecular environment available for this source do not allow us 
to make conclusions about the nature of the 9.9-GHz maser. As 
previously mentioned in section lJTl there is an Si O emission (pre- 
sumably tracing shocks) detected in this source jPe Buizer et al.l 
1 200ft) at a position which does not allow us to completely rule 
out the association with the maser. On the other hand, the 9.9- 
GHz maser is located near the protrusion on the southern side of 
the Hit region (see Fig.[2ji) which may reflect the presence of the 
same type of interaction between the ionised and neutral material 
as described above. Although the remaining source G343.12— 0.06 
is a clear cas e where the 9.9-GHz m aser is associated with an out- 
flow (Fig.[3h: l\tonkov et alj200 ?). these considerations show that 
the outflow scenario may not be the only option. 

This idea should equally apply to other class I masers, includ- 
ing 36 and 44-GHz ones, which are more ubiquitous than the 9.9- 
GHz masers. Therefore, future interferometric surveys of class I 
masers, e.g. a follow up of the MMB survey (briefly mentioned in 
ICasweU et aT]|2010l) . have a great potential to test this hypothesis. 
It is worth mentioning that lKogan & Slvshl ( 1 19981) observed W33- 
Met at 44 GHz. However, phase referencing was not used in their 
experiment making the search for possible associations difficult due 
to a lack of absolute positions. 

4.4 Evolutionary stages 

The question whether different masers trace distinct evolu- 
tionary stages of high-mass star for mation has been investi- 
gate d for m ore than two decades (e.g., ' Forster & Caswelilll989l; 
■ Codella & M oscadelli 2000; Ellingsen 2005l l2006l ; IVoronkov et al] 
120061 ; IBreen et alj 120 1 Q) . The majority of class II methanol 
masers at 6.7 G Hz are associated with millimetre and submil- 
limetre sources jWalsh et alj |2003| ; iHill et alj |2005|) . while only 
a small fraction shows association w ith detectable Hll regions 
jPhillips et al]|l998l : IWalsh et alj|l998h . This leaves no doubt that 
they trace a very early stage of the (proto)stellar evolution preced- 
ing the formatio n of an Hll region, and fade out shortly after an 
Hll region forms. ISreen et alj ( l201(]|) find that the other widespread 
class II methanol masers at 12.2 GHz most likely lag behind the 
6.7-GH z masers. The latter se em to largely overlap with the H2O 
masers jSzvmczak et al.ll2005l) whereas the OH masers are usually 
considered to be a signature of a somewhat more evolved stage, 
since they show a greater overlap with the ultrac ompact Hll regions 
than is typical for class II methanol masers jForster & Caswell 
1989; Caswell 1997). 

However, it is still poorly understood as to where the class I 
methanol masers fit into this picture. The main obstacle is the lack 
of unbiased surveys of the class I masers, which cover a substan- 
tial part of the Galaxy. As discussed above, the majority of cur- 
rently known class I methanol masers have been found in sur- 
veys targeted towards class II methanol masers. Therefore, the 
sample of currently known class I masers is strongly biased to- 
wards the evolutionary stage traced by the 6.7-GHz masers. With- 
out a widespread and bright low-frequency class I methanol maser 
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transition (similar to the class II maser transition at 6.7 GHz or 
mainline OH) a bl ind survey is very time consuming. Therefore, 
iPratap et al ] (|2008) adopted an interesting strategy to cover sev- 
eral known molecular clouds in their, otherwise untargeted, search 
for class I methanol masers at 36 and 44 GHz. Com parison with 
the pr esence of Hll regions and infrared sources, led IPratap et al.l 
1 I2OO8I) to suggest that the 36-GHz masers may appear earlier than 
the 44-GHz masers during the evolution of a star-forming region. 
However, these data do not allow interrelations between evolution- 
ary stages where t he class I and cla ss II methanol masers are present 
to be established. lEllingsenlJioO^ investigated this question using 
the infrared properties of a subsample of methanol masers asso- 
ciated with GLIMPSE (Galactic Legacy Infrared Mid-Plane Sur- 
vey Extraordinaire) catalogue point sources. It was found that the 
sources with a class I maser seem to have redder GLIMPSE colours 
than those without it. With the caveat that the small number statis- 
tics made any firm conclusion impossible, EUingsen (20Qi) sug- 
gested that the class I methanol masers may signpost an earlier 
stage of high-mass star formation than the class II masers. 

These and other considerations led lEllingsen etall j2007l) 
to propose a qualitative evolutionary scheme for the different 
maser trans i tions, which has subsequently been further refined by 
iBreen et al.l ( 1201 Oh in their Figure 6. In this sequence the class I 
masers represent the first signatures of high-mass star-formation 
and fade out within approximately 2x10* years after the onset 
of the 6.7-GHz maser emission . Base d on the detailed study of 
G343. 12-0.06, IVoronkov etal] ( |2006|) pointed out that the sug- 
gestion that class I methanol masers precede class II is inconsis- 
tent with the presence of OH masers and the lack of the 6.7-GHz 
maser emission in this source. OH maser emission in the absence 
of 6.7-GHz masers is more consistent with an evolved stage of star- 
formation (see e.g. , Caswell 1997) . 

It is important to keep in mind the major ass umptions which 
underl ie the evolutiona ry tim eline put forward by EUingsen et al. 
( I2OO7I) and Breen et al] j2010l) . In particular, that each major maser 
species arises only once during the evolution of a particular star 
formation region, and that all the maser species are associated with 
a single astrophysical object. For regions where the positions of 
the different maser species are known to an accuracy of an arc- 
second or better this latter assumption can be proven or disproved 
for most transitions; however, the larger spatial (and angular) scale 
over which class I methanol masers are often observed makes it 
muc h less clear whether the emission is due to a single object (see , 
e.g., 'Kogan & Slvsh"l998l; iKurtz et alj2"004l ; IVoronkov et alj2006l : 
ICyganowski et al. 2009) . The degree of ov erlap of the class I and 
class II methanol masers on the timeline o f lEllingsen et"ai] j2007t) 
and lBreenet"ZI ( l20ld ) is based on the rate of detection of class I 
metha nol masers towards class II masers observed by lEllingser] 
feoosi). However, the idea that class I methanol masers precede 
the class II was based purely on the circumstantial evidence that 
the former are associated with outflows, which are generally con- 
sidered to be signposts of young objects, and the suggestion of 
[EUingsen (2006) that the class I masers may be associated with 
redder GLIMPSE point sources. As discussed earlier, all the 9.9- 
GHz masers found in this study are located near Hll regions (see 
section |431 >. all but W33-Met have OH masers detected (see sec- 
tion l3.1l no sensitive OH data seem to be available for W33-Met), 
and one source, W33-Met, does not have a 6.7-GHz maser de- 
tected (see section l42l (. Despite the small number of sources, these 
similarities encourage us to revisit the question of where class I 
methanol masers fit into a maser-based evolutionary sequence. 

Contrary to the sequence presented by iBreen et^ i2010t) it 



is clear that class I masers do overlap in time with OH masers. 
The 9.9-GHz maser detections reported in this paper reinforce 
the conclusion, although the class I masers in o ther transitions 
were already known for all of these sources (e.g ., iBachiller et al.l 
I1990I ; ISlvshetalJlT994l : iLiechti & Wilsorilll996n . Despite the bi- 
ased target selection in the majority of class I maser surveys, a 
few sources such as W 33-Met (see section lSTt and G343.12— 0.06 
dVoronkov eTZIb OOe). are known to have class I masers without a 
class II counterpart at 6.7-GHz (see also section |42l l. This suggests 
that it is more appropriate to place the class I methanol masers as 
partly overlapping, but largely post-dating the evolutionary phase 
associated with class II methanol masers. As a part of the follow up 
program b ased on the MMB untargeted survey briefly mentioned 
by Ca swell et alj JioTQ) . class I methanol masers at 44 GHz were 
searched for towards a number of OH masers which have no de- 
tected 6.7-GHz masers. The full results of this work will be re- 
ported elsewhere. However, the preliminary analysis indicate a de- 
tection rate greater than 50%. This supports the hypothesis on the 
lifetime of class I masers given above and moreover suggests that 
the result is not unique for the 9.9-GHz masers and also applies to 
the widespread class I methanol masers (e.g. that at 44 GHz). 

The class I masers are formed in moderately dense gas 
separated from stron g infrared sources (for discussion see, e.g., 
IVoronkov et"^l2006h . These masers have radial velocities close 
to that of a parent molecular cloud which suggests, together 
with morphology, an association with the interface region. Weak 
shocks interacting with the molecular cloud dramatically change 
the chemical composition of the medium by releasing methanol 
and other complex molecules from the dust grain mantles into the 
gas phase. This increase of methanol abundance is observed in 
the sh ock-processed regions related to outflows (e.g. lGibb & Davii 
Il998l) . However, the same effect is also expected if other mech- 
anisms, not just outflows, drive shocks into the parent molecular 
cloud, provided the shocks are weak enough and do not dissoci- 
ate methanol molecules. In particular, as discussed in section 143] 
class I methanol masers could be associated with expanding Hll 
regions as well as with outflows. The presence of an Hll region im- 
plies that the star formation in a particular site is quite evolved. On 
the other hand, an ample amount of molecular material interacting 
with an outflow or expanding Hll region is required to produce a 
detectable maser. Therefore, there could be a bias towards deeply 
embedded regions of star formation. The amount of quiescent ma- 
terial remaining in a star-forming region showing maser activity de- 
pends not only on the age, but also on the initial conditions which 
may vary from one site to another. 

It it important to recall that high-mass star formation usually 
occurs in a crowded environment where sequential or triggered star 
formation may take place. G19.61— 0.23 and W33-Met are good 
examples of clustered star formation. For example. IHo et al.l l fl986l) 
suggested that the interaction of ionised gas from an earlier epoch 
of star formation may have triggered the formation of a second 
generation stars in W33-Met through the radiation-driven implo- 
sion mechanism. Therefore, in some cases it may be impossible to 
attribute a class I maser spot to a particular YSO, especially tak- 
ing into account that the angular offset of this maser may exceed 
the separation between YSOs in the cluster. In addition, the indi- 
vidual constituents of the cluster may have significantly different 
ages. For example, one YSO in the cluster may have an associ- 
ated 6.7-GHz maser, while another could be associated with OH 
masers and a well developed Hll region. The fact that the class I 
methanol masers can be associated with both outflows and expand- 
ing Hll regions, which presumably appear at different evolutionary 
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stages, further complicates the situation. At this stage the number 
of sufficiently simple sources studied in detail is quite limited hin- 
dering the statistical analysis. The future interferometric surveys of 
the 44- and 36-GHz methanol masers (widespread class I maser 
transitions) will be vital for refining the evolutionary sequence for 
masers. 



4.5 Association with tlie EGOs 

IChen et alj ( l2009h fo und that approximately two-thirds of the 
EGOs from the list o f ICvganowski et alj feOOsh . which have been 
observed at either 95 or 44 GHz, have a class I methanol maser 
detected within 1'. Based on this result, they made a comment that 
EGOs could be good targets for future class I m aser searches. Such 
a surv ey of the 44-GHz masers carried out by ICyganowski et"ail 
indeed had a high detection rate (90%, 17 detections out 
of 19 EGOs observed). However, two out of four 9.9-GHz masers 
reported in this study are not associated with an EGO (see sec- 
tion I3.lt . The physical nature of EGOs requires a better under- 
standing dCyganowski et alj2009l) . One may speculate that the out- 
flows are more likely to produce an EGO when they interact with 
the parent molecular cloud, than if shocks are caused by an ex- 
panding Hll region. It is worth mentioning, that the survey of 
ICvganowskietalJ ( l2009l) discovered an interesting class I methanol 
maser G49.27— 0.34. No 6.7-GH z mas er has been detect ed in this 
source dCyganowski et al.ll2009l) . while lMehrin"gerl h994h detected 
smooth radio continuum from presumably a well evolved Hll re- 
gion (seen at 20cm only). 



5 CONCLUSIONS 

(i) Two new class I methanol masers at 9.9-GHz were found 
in addition to two other sources already reported in the literature. 
The absolute positions with arcsecond accuracy are summarised 
in Table |3] for all four 9.9-GHz masers known to date. Based on 
the trial observations, we also suspect that another 9.9-GHz maser 
may exist in G305.21+0.21 at a flux density below the detection 
threshold of our survey. 

(ii) We suggest that some class I methanol masers may be asso- 
ciated with shocks driven into molecular cloud by an expanding Hll 
region. This is an alternative scenario to the association with out- 
flows which has been proved unambiguo usly in a n umber of other 
cases (e.g.,, Kurtz et al. 2004 ; Voronkov et'SfeOOel) . This new sce- 
nario applies also to the class I methanol masers in other transitions, 
e.g. to the widespread masers at 36 and 44 GHz. It does not appear 
to be an exclusive mechanism to generate the 9.9-GHz masers, al- 
though could be responsible for 3 out of 4 known 9.9-GHz masers. 

(iii) The evolutionary stage with the class I maser activity is 
likely to outlast the stage when the 6.7-GHz methanol masers are 
present and overlap significantly in time with the stage when the 
OH masers are active. We expect that the class I masers in one 
of the widespread maser transitions (e.g. at 44 GHz) will be de- 
tected towards a significant number of OH maser sources which do 
not have a class II methanol maser at 6.7 GHz. Such sources were 
rarely observed in the class I maser surveys currently available in 
the literature. 

(iv) A 9.9-GHz methanol maser was found in G331.13— 0.24 
which also hosts one of the few known periodically variable class II 
methanol masers at 6.7 GHz. This rare combination makes this 
source an attractive target for a long-term monitoring program. 



which has a high potential to shed light on the periodic variabil- 
ity of masers. 
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